1 



Unquenched domain wall quarks with TSMB 

qq+q Collaboration 

F. Farchioni a , C. Gebert b , I. Montvay b *, L. Scorzato b 



(N 
O 
O 
(N 

< 

00 
(N 



> 

O 
oo 
O 
<N 
O 

in 



X 



a Institut fiir Theoretische Physik, Universitat Munster, Wilhelm-Klemm-Str. 9, D-48149 Munster, 
Germany 

b Deutsches Elektronen-Synchrotron, DESY, Notkestr. 85, D-22603 Hamburg, Germany 

The numerical simulation of domain wall quarks with the two-step multi-boson (TSMB) algorithm is investi- 
gated. Tests are performed on a 8 3 • 4 lattice with Nf — 2 quark flavours. 



1. INTRODUCTION 

Lattice actions with improved chiral properties 
offer the perspective of QCD simulations with a 
better control of chiral symmetry at non-zero lat- 
tice spacings. A prototype is the domain wall 
fermion action [|l}-|| defined on a five dimensional 
hypercubic lattice. Following Shamir's formula- 
tion ||, the light chiral fermion modes are lo- 
cated on two boundaries of the fifth dimension. 
The price of the chiral symmetry at non-zero lat- 
tice spacing is the extra dimension enlarging the 
number of degrees of freedom. From a technical 
point of view this means that the extensions of 
the fermion matrix are larger and therefore the 
numerical simulations are slower. 

It is an interesting question how much the com- 
putation speed decreases compared to, say, Wil- 
son fermions. Since the good chiral properties 
of domain wall fermions develop only sufficiently 
close to the continuum limit, the performance 
studies have to be finally performed on large lat- 
tices. A comparison can only be conclusive if the 
autocorrelations of important physical quantities 
are also determined - a difficult task requiring a 
substantial amount of computer time. 

A first step towards the effective simulation of 
domain wall quarks is to identify possible sim- 
ulation algorithms which have the potential of 
being applicable in simulations with small quark 



masses and in large physical volumes. The two- 
step multi-boson algorithm (TSMB) g-0 has 
been succesfully tested in this respect in case of 
the Wilson quark action |^|| . The application of 
TSMB for domain wall quarks has been recently 
considered in |lO[ |. Here we report on some fur- 
ther tests along these lines. 

2. ACTION AND ALGORITHM 

The lattice action for domain wall quarks can 
be introduced as 

S = S G [U] + S F [V, , U] + S PV [&,$, U] . (1) 

Here Sq denotes the pure gauge-field part de- 
pending on the gauge field U , Sf is the fermionic 
part with the Grassmannian quark fields "J, ^ and 
Spy depends on the bosonic Pauli-Villars field $ 
which subtracts the heavy fermion modes - as in- 
troduced in p[ . 

The fermion action Sf is defined by 



S F = 



x.s; x' ,s' 



V(x',s')D F (x',s'; x,s)V(x,s) . (2) 



The four-dimensional space-time coordinates are 
denoted by x, x' and the fifth coordinates are 
s, s' . The domain wall fermion matrix Dp is con- 
structed from the standard four-dimensional Wil- 
son fermion matrix 
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D(x', x) = b X ' X (4 - am a ) 
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+ 4'+A,x(1-7m)^] 



(3) 



The notations are standard: a is the (four- 
dimensional) lattice spacing and fl denotes the 
unit vector in direction fi. The bare mass —mo is 
chosen to be negative and should be tuned prop- 
erly for producing the light boundary fermion 
state. The non-zero matrix elements of the do- 
main wall fermion matrix are: 

{D F ) S . S =a + D , (D F ) S!S+1 = -aP L , 
(D F ) s+1;a = ~aP R , (D F )x >Ns = am f P R , 
(D f )n„i = am f P L . (4) 

Here to/ denotes the bare fermion mass of the 
light boundary fermion, N s is the lattice exten- 
sion in the fifth dimension, Pr = 5(1 + 75), Pl = 
4(1 — 75) are chiral projectors and a = a/a s de- 
termines the lattice spacing in the fifth dimension 
a s relative to a. 

The TSMB algorithm is based on the her- 
mitean domain wall fermion matrix 

D F = 75 i? 5 D F (5) 

where the reflection in the fifth dimension is de- 
fined as (R 5 ) s '.s = 5 s ', Ne+ i- s (1 < s,s' < N s ). 
The quark determinant with Nf flavours is real- 
ized in TSMB by 

N f 1 

.(6) 



(2), 



det (£>*■) 

det P£>(D F ) detP^'(Dp) 
The polynomials P^ and P^ satisfy 
lim pW(x)Pg>( x )= x -"f/z , x £ [e, A] 



(7) 



The interval [e, A] covers the spectrum of Dp on a 
typical gauge configuration. The first polynomial 
PW of order n\ is a crude approximation and is 
realized by the multi-boson representation |12"| . 
The second polynomial P^ of order ^> «i 
is a better approximation. It is taken into ac- 
count in the updates by a global noisy Metropo- 
lis correction step. Since for fixed 712 (and outside 
the interval [e, A]) the approximation p( 1 )p( 2 ) is 
not exact, a final correction step is performed by 
reweighting the gauge configurations which are 
considered for the evaluation of the expectation 
values. (For mode details on TSMB see H-@.) 



3. SIMULATION TESTS 

The implementation of TSMB for domain wall 
quarks is straightforward. Since domain wall 
fermions have an extra index labeling the fifth 
coordinate, a potential problem is the storage of 
the ri\ multi-boson fields in computer memory. 
This problem can be easily solved because (for 
domain wall fermions and in general) one can or- 
ganize the gauge field update in such a way that 
the dependence on the multi-boson fields is col- 
lected in a few auxiliary 3 ® 3 matrix fields which 
can be easily stored in memory. The multi-boson 
fields themselves can be kept on disk and have to 
be read before and written back to disk after a 
complete boson field update. The duration of the 
input-output is negligible compared to the time 
of the update. Organized in this way, TSMB has 
a rather low storage requirement. 

Wave function profiles in the fifth dimension 
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Figure 1. Wave function profiles for eigenstates 
of the hermitean domain wall fermion matrix. 

We performed test runs for two degenerate 
quark flavours Nf = 2 on 8 3 • 4 lattices in the 
vicinity of the N t = 4 thermodynamic crossover. 
The parameter sets have been chosen from the 
points in parameter space which were investi- 
gated in |1^| . Typical parameters were: M = 
cztoo = 1.9, M q ee anrif = 0.1, a = 1.0, 0.5 and 
5.20 < 13 < 5.45. (The Pauli-Villars mass was 
Mi = 1.) After equilibrating the gauge configu- 
rations several features of the domain wall quarks 
were investigated. For typical wave function pro- 
files in the fifth dimension with N s — 8 see fig- 
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ure [l]. The curves peaking at the walls corre- 
spond to eigenvalues with smallest absolute value, 
those concentrated in the middle to eigenvalues 
with largest absolute value. Spectral properties 
of fermion matrices are shortly discussed in the 
next section. 

4. SPECTRAL PROPERTIES 

As is well known from quenched studies (see 
for instance fll4| ), the good chiral properties of 
domain wall fermions are realized only if the 
four-dimensional fermion matrix D used for the 
construction of the domain-wall fermion matrix 
Df does not have very small eigenvalues. The 
hermitean four-dimensional fermion matrix D = 
75 D should have a "gap" near zero in its spec- 
trum. 



Spectral flow with the Wilson-Dirac fermion matrix 



Spectral flow with the domain wall fermion matrix 
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Figure 2. Eigenvalue flow for the hermitean 
Wilson-Dirac matrix. 

In our unquenched test runs such a gap does 
not appear (see figure ^J). The lattice spacing is 
obviously still too large. The spectrum flow of 
Dp on the smae gauge configuration is shown by 
figure ||. 

The conclusion of our first tests is that the ap- 
plication of the TSMB algorithm for numerical 
simulations of domain wall quarks is straightfor- 
ward. A comparison of the computation speed 
compared to, say, Wilson quarks requires an anal- 



11.2 



0.1 



>• 



-0.2 







-o.i 



-0.05 



0.05 



0.1 



Figure 3. Eigenvalue flow for the hermitean 
domain wall fermion matrix. 

ysis including the measurement of autocorrela- 
tions. Since the good chiral properties of domain 
wall fermions develop only sufficiently close to the 
continuum limit, the performance studies have to 
be carried out on large lattices. 
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